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We derive the threshold-resummed total cross section for heavy quark production in hadronic 
collisions accurate to next-to-next-to-leading logarithm, employing recent advances on soft anoma- 
lous dimension matrices for massive pair production in the relevant kinematic limit. We also derive 
the relation between heavy quark threshold resummations for fixed pair kinematics and the in- 
clusive cross section. As a check of our results, we have verified that they reproduce all poles of 
the color-averaged qq — > ti amplitudes at two loops, noting that the latter are insensitive to the 
color-antisymmetric terms of the soft anomalous dimension. 

I. INTRODUCTION 

Heavy quark production has a long and continuing history of phenomenological interest. Top-pair production is 
among the most important standard model processes in the context of searches for new physics. It is also among the 
most challenging in computation. Indeed over 20 years passed between the landmark first NLO calculations of heavy 
quark production [l[ and the derivation of analytic expressions for the inclusive cross section. 

As in many hard hadronic processes, higher order perturbative calculations contain logarithmic enhancements, 
associated with the approach to partonic threshold, that is, configurations where the initial-state partons have just 
enough energy to produce the observed final state. Threshold resummation Q organizes such logarithms, in a manner 
we review below. The current state of the art for heavy quark production is next-to-leading log (NLL) resummation 
matched onto exact next-to-leading order (NLO) results. In this paper, we study the extension of this formalism to 
next-to-next-to-leading log (NNLL) resummation and beyond in the production of heavy quark-antiquark pairs at 
hadron colliders. In this context, we derive the NNLL inclusive cross section, generalizing the NLL results of [H, 
starting from the resummation formalism at fixed parton kinematics described in [f| . The numerical impact of our 
present work will be detailed elsewhere. 

While our considerations are reasonably general and can be applied to a larger class of processes, we choose to 
present our results and their derivation for heavy quark production because recent advances [H 0, HI make it possible 
for us to determine the essentially two-loop inputs necessary for explicit NNLL resummation in this case. These inputs 
are the so-called soft anomalous dimension matrices for heavy quark pair production 5] which we exhibit below in the 
relevant kinematic configuration for the total cross section. Drawing from the formalism of Ref. j|, we will perform 
a single one-loop calculation, corresponding to a boundary condition in the evolution of soft gluon emission, which is 
necessary for the complete NNLL result. As noted at the level of NLL, for the inclusive cross section, resummation 
can be carried out separately for pair production in the s-channel color singlet and octet states, without color mixing. 
We will find the same structure at NNLL, and present our final results in a form that follows the NLL formalism of 
Ref. [|. 

We begin with a review of threshold resummation for semi- inclusive cross sections at fixed kinematics for the partonic 
scattering process, in the formalism developed in Ref. and applicable in principle to any logarithmic approximation. 
This formalism relies on the factorization of color-diagonal "jet" functions associated with the external energetic and/or 
massive partons that take part in the scattering at short distances and a "soft function" describing the exchange of 
low-energy quanta between these particles. We identify in particular, a scheme to resolve the ambiguity between the 
jet and soft functions, based on the singularities of form factors in dimensional regularization. We go on in Section Hill 
to show how an expression for the threshold-resummed inclusive cross section may be derived from the resummation 
at fixed kinematics. With this result in hand, we determine two-loop soft anomalous dimension matrices in Sec. IIV1 
which we need to determine the soft functions at NNLL. We describe a very nontrivial check of these results, which 
fully reproduce the two-loop pole structure of heavy quark pair production in light quark annihilation. In Sec. fVj we 
assemble the remaining ingredients in the NNLL resummation, including the one-loop boundary condition mentioned 
above. We conclude with a summary and a few comments on prospects for future work. 
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II. THRESHOLD RESUMMATION AT FIXED KINEMATICS 



In this section, we review the threshold resummation formalism of Ref. [5j], which is adapted to semi- inclusive 
reactions characterized by fixed partonic scattering kinematics, as in for example, 

A(P1) + /a(P2) - fa(Pa) + fbipb) , (1) 

where fi(pi) denotes a parton of flavor fi and momentum pi. We have shown a 2 — ► 2 process, but final states with 
more than two particles are also possible, so long as all invariants pi ■ pj are large. The formalism we sketch in this 
section applies to processes involving light quarks and gluons, and also to the production of heavy quarks. In the 
latter case, we can also study the inclusive cross section, for which threshold resummation has been developed from a 
related point of view Q • In Section IIIII we will derive resummed inclusive cross sections for heavy quark production 
from their semi-inclusive forms. 



A. Factorization near partonic threshold 

Our starting point for the resummation of observables involving initial and/or final state hadrons is the formalism 
of Ref. To be specific, we restrict our discussion to the 2^2 processes of Eq. (p}, although many of our 
considerations can be directly generalized. For the production of a pair of particles with mass m, the kinematics can 
be described by the invariant mass M and rapidity y of the partonic final state and the pair center-of-mass rapidity 
difference r). Assuming that m 3> Aqce>, this cross section can be written in standard factorized form as 

T, / ri dtJ h 1 h 2 ^QQ f 1 , f dx a dx b 2 2 



x 8 I z '■ — I 5 ( y — - In X " 



M 2 m 2 



f 



where we have normalized the cross section so that all quantities are dimensionless. The purpose of threshold 
resummation is to organize plus distributions in the variable 



M 2 



X a Xb X a XfoS 



(3) 



with x a and x b the usual fractional momenta. Partonic threshold is defined as the limit z — * I, at which the incoming 
partons provide just enough energy to produce the observed final state. The mismatch between real gluon emission 
and virtual corrections gives rise to singular distributions at z = I. These distributions appear in the nth order 
expansion of the perturbative function lo^j^qq up to the level of a™[ln 2n_1 (l — z)/(l — z)] + . 

In Ref. it was observed that as z — > 1, partonic cross sections can be factorized into a set of universal factors 
associated with the incoming and outgoing partons of the underlying process, along with process-dependent factors 
that describe the coherent interactions of those partons, at short and long distances. The resummed dependence in 
I — z is conveniently generated by taking Mellin moments with respect to z, schematically, 

a{N) = / dzz N - 1 a(z) 



o 



l 

dze^ l - N - 1){1 - z) a{z)+0{l/N). (4) 



o 

For example, in the inclusive Drell-Yan process, the corresponding kinematical variable is z = Q 2 /s, where s is the 
partonic cm. energy [3| . For resummation of the total inclusive cross section of heavy quark pair production at hadron 
colliders the corresponding variable is z = Am 2 / s where m is the mass of the top quark. Another example relevant 
for this paper is the invariant mass M 2 ^ = {p t + pi) 2 distribution of a top quark-antiquark pair; the relevant partonic 
variable is z = M 2 ^/ s. In any case we assume that the partonic variable is defined such that threshold kinematics is 
attained in the limit z — * 1. In moment space this corresponds to the limit N — * oo. The analysis of Ref. [5( exploits 
factorization near threshold, according to which the cross section can be written as a convolution in an appropriate 



3 



momentum component of the soft radiation associated with a set of functions 0. In threshold resummation for 
hadronic collisions, this component is the energy, E* , of each final-state particle in the center-of-mass frame of the 
hard collision. That is, for any threshold resummation at hadronic collisions, we can identify 



E 



particles 



2EI 



(5) 



where the partonic variable s = x a XbS equals M 2 at threshold, with M the invariant mass of the observed pair of heavy 
particles. The cross section then factorizes into simple products in the corresponding moment space. Dependence on 
the moment variable enters only through the transform, and is therefore always in the form N/M, up to corrections 
that decrease as powers of N. 

As a result of this analysis, the partonic cross section takes a factorized form in moment space, which we can 
represent as 



ujp 



M m 

N,fj,— ,— ,a s (/i 2 ) ) = Ji(iV,a s (/i 2 ))... J k (N,M/n,m/fi,a s (Li 2 )) 



x Tr 



LI LL 



N 2 LL 2 M 



M 2 



0(1/N) , (6) 



where the label P refers to a particular partonic process, for example qq — ► ti, with q a light flavor. The Mellin 
moment N is conjugate to the kincmatical variable z. As shown, the various functions appearing in Eq. ([B]) depend 
on other kinematical variables and masses as well as the factorization and renormalization scales. These functions 
depend on the specific process. Below, we will give them more explicitly in the specific examples considered here. We 
will refer to the factors Ji appearing in Eq. ^ as the jet functions for the underlying process. They are color diagonal 
functions that describe the factorized dynamics of initial and/or final state hard partons, whether massive or massless, 
and as such are independent of the details of the hard subprocess. Jet functions for initial-state partons absorb the 
collinear subtractions necessary to define the hard scattering function u in Eq. ([S]), so that they are infrared safe. Jet 
functions for final-state partons are automatically infrared safe for the differential and inclusive cross sections that 
we discuss here. The formalism can be extended as well to a variety of jet observables and to single-hadron cross 
sections. The number k of such functions in Eq. ([6]) corresponds to the number of hard colored partons in the process 
being considered. 

The functions H and S appearing in Eq. ^ are known as hard and soft functions, respectively. They are both 
matrices in the space of tensors that describe the exchange of color at short distances |5j. Examples for quark- 
antiquark scattering are color singlet or octet in the s- or i-channel. We will denote these tensors in boldface, and 
their product is traced over the combinations of color tensors in the amplitude and its complex conjugate. In the limit 
N — > oo the hard function H is free of logarithmic dependence on N; it is obtained from a dedicated, process-specific 
calculation. 



B. Moment-dependence and the soft anomalous dimension matrix 

The soft function S contains terms due to wide-angle soft emissions and thus contributes a single power of hi(N) 
per loop. It is also process dependent, and in the general case is dependent on the four-velocities {/3i} of the partons 
that take part in the hard scattering. For processes involving four or more colored hard partons it is a matrix in 
the space of color tensors. Assuming fixed-angle scattering, the soft function depends on the scalar products of these 
velocities, in addition to a single overall scale, which we will take to be M, the invariant mass of the pair for the 
case of heavy quark production. For a massive quark of velocity (3 q , we shall set f3 2 = m 2 /M 2 , and for most of this 
discussion, treat this ratio as a number of order unity. 

As noted above, all iV-dependence is of the form N/M. As a result, in the dimensionless soft function, iV-dependence 
appears only in the combination M/{Nli). In Ref. [H, it was shown that the TV-dependence of the soft function 
S(N, . . . ) entering the cross section Eq. © can be made explicit in terms of a "soft anomalous dimension matrix", 
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Tg. Making the natural choice, fi — M, we have 

- R, . R: nln 2 \\ — Vwn) - I 

M 2 



= V exp ■ 



^4 fa. factf)) 

M/N V 



xS(l,ft -ft,a s (M 2 /N 2 )) 




o 



M ^r s (A-ft,a s (// 2 ))} 

M/N f 1 J 
1 T N-1 _ 1 . 

<*c l _ x r+ (A ■ /J,-, a s ((1 - x) 2 M 2 )) 



xS(l,ft -ft,a s (M 2 /N 2 )) 




N-l _ i 1 

— — T S (A: • ft,a s ((1 - x) 2 M 2 )) I , (7) 

where the second expression is accurate to next-to-next-to leading logarithms (i.e. terms ~ a™ In" N in the cross 
section) for N = iVe 7E , with je the Euler constant. Throughout this paper a s = a s (n 2 ) is the standard MS coupling 
evolving with Nl light flavors. Decoupling of the heavy flavor will simplify our results significantly. The relation 
between the bare a b s and renormalizcd couplings reads 



,b 



ft = a s (n 2 



i _/3o^) + ' 
4e 7r 



(8) 



where S e = (47r) e exp(— £7e) and /3 = (11/3)Ca — (4/3)TjpNl. The color factors in an SU(iV)-gauge theory are 
C A = N,C F = (N 2 - 1)/(2N) and T F = 1/2. 

The structure of Eq. (0 follows from the renormalization group equation satisfied by the soft function 
S{N 2 y 2 / M 2 , . . . ), where Tg plays the role of a matrix of anomalous dimensions The function S(l, . . . ) plays 
the role of a boundary condition, which is chosen to be the soft function at unit N, that is, with unit weight. In 
general, this factor contributes a single ln(iV) starting from two loops, which is due, however, entirely to the presence 
of N in the scale of the running coupling in its one-loop expression. To determine this contribution one need only 
calculate the soft function in Eq. ([6]) through one loop. 

At N — 1, the computation of the soft function is given by a total eikonal cross section, subtracted for eikonal 
jet functions to eliminate collinear enhancements [f|. In the formalism of Ref. virtual corrections are pure 
counterterms, because the corresponding eikonal diagrams are scaleless and vanish in dimensional regularization. In 
the full soft function, however, the hard scale sets a maximum total energy for the soft function at N = 1, and 
the corresponding integrals are not scaleless. Their infrared poles are cancelled by the virtual diagrams, but finite 
corrections may remain. 

In summary, the soft function S at N = 1 takes the form 

S (1, ft • ft, a s (M 2 /N 2 )) = S<°> + ^ ( M l' N2) S« (1, ft • ft) + . . . , (9) 

where S' ' is a constant diagonal matrix independent of the coupling and S^) (l,ft • ft) is free of dependence on N, 
but can depend on the eikonal velocities that define the soft function. Explicit expressions for relevant to heavy 
quark production can be found in [lol ]. We will give the one- loop correction below, after specifying a scheme that 
defines the soft function unambiguously. At this stage, we note that to compute the soft function fully at next-to- 
next to leading logarithm it is necessary to compute the two-loop anomalous dimension matrix and the one-loop soft 
function. 



C. The form factor scheme 



The soft function is not unique, but is ambiguous at the level of single logarithmic contributions that can be 
absorbed into the jet functions. These ambiguities, must be proportional to the unit matrix in the color exchange 
space (since the jet functions are diagonal in color). To resolve this ambiguity one has to specify a prescription for 
the definition of the anomalous dimension matrix r<j, which we discuss next. 
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A fundamental observation of Ref. @ is that the matrix T$ appearing in Eq. {7J can be extracted from the 
corresponding (virtual) amplitude for the process under consideration. To review how this can be done, we first 
observe that any on-shell, renormalized scattering amplitude at fixed angles can be factorized as follows [111 ] : 

M/(e, . . . ) = h(e, ...).. .j s (e, ...) ■ hj{. ..)■ s u (e, ...). (10) 

Here I, J are indices that label color exchange tensors; in particular, they indicate that the amplitude M can be 
thought of as a vector in the space of color representations In order to emphasize the similarity between the 

objects appearing in Eqs. ([6]) and ([T0|) we have used the same letters to denote jet, soft and hard functions. It 
should be stressed, however, that these are not the same objects. In particular, the moment N does not appear in 
an amplitude. To distinguish clearly between the objects appearing in the cross section and in the amplitude we use 
lower case letters, and explicitly show the dependence on the infrared regulator e where d = 4 — 2e. 

For both massive and massless external partons the amplitude soft function sjj(e, . . . ) appearing in Eq. (|10j) is fully 
determined by its (matrix) anomalous dimension Tjj. It can be computed order by order in perturbation theory, as a 
series in the couplingo; s (/i 2 ). The properties of Tjj have been studied extensively in the massless [H, E3, EE E3] an d 
in the massive [g, 0, Q cases. To unambiguously fix the soft function in the massless case, a natural and convenient 
scheme was proposed in Ref. [13] • There, the jet functions for each external parton were identified with the square 
root of the massless on-shell space-like form factor for the corresponding parton. We will assume this prescription, 
which we call the form factor scheme by default from now on. In this scheme, the jet functions for massless particles 
are series in the coupling a s (^L 2 ) with coefficients that are e-dependent numbers. A natural extension [181 ] in the 
massive case is to identify the jet function with the small-mass limit of the corresponding massive space-like QCD 
form factor. Therefore, in the massive case the jet function contains also powers of ln(/i 2 /m 2 ), where m is the pole 
mass of the heavy quark. 

In the form factor scheme, we can derive an explicit expression for s(e, . . . ) in terms of the matrix soft anomalous 
dimension, 



s(e, . . . ) = Pexp | - jf 1 J^- T S (a s [(1 - x) 2 Q 2 } ) j , 



(11) 



where T$ = (a/ir)Tg + (a/Tr) 2 T^ + Ofa?), and a stands for either a s (/i 2 ) or a s . The coupling a s (k 2 ) is the d- 
dimensional strong coupling constant [HI, |2C| . known through NNLO [2l[ . It is a function of the usual four dimensional 
coupling a s (fj, 2 ) and the dimensional regulator e (the explicit relation we use here can be found in Ref. [IH [60]). The 
result of Eq. (llip depends on the vanishing of the running coupling a s (/x 2 ) at [i = for e < 0, that is, in more than 
four dimensions. 

Note the similarity between Eq. (jlip and ([7]) : the amplitude term can be directly obtained from one of the exponents 
in Eq. |(7J) by simply ignoring the term with a;^ -1 and replacing the four dimensional coupling with the d-dimensional 
one. Such duality is not accidental; physically both exponents can be thought of as two different rcgularizations of 
the soft limit with regulators, respectively, e and ln(iV). This relation has been explored and detailed, for example, 
in Ref. QJ]. 

Without affecting the value of the soft anomalous dimension matrix, one has the freedom to add finite e-terms in the 
soft function which amounts to a re-definition of the hard function h. In Eq. (Illj) we choose a minimal, MS-inspired 
scheme, where only e poles are kept in the soft function. With this scheme defining the separation between the soft 
and hard functions, the explicit relation between the soft function and the anomalous dimension matrix through two 
loops reads: 

To summarize our discussion up to here, we have shown that the soft function in Eq. ([6]) can be fully specified 
by (|7l9p in terms of the anomalous dimension matrix Tg which in turn is derived solely from the knowledge of the 
purely virtual corrections to the same process. Ambiguities in Ts are fixed by choosing a prescription at the level 
of the amplitude and we work with the form factor prescription for both massless and massive partons. Once the 
soft function in Eq. ([6]) has been fixed, in order to perform NNLL resummation in observables, one has to determine 
the number and the form of the various jet functions related to that observable in a manner consistent with the 
prescription implicit in the definition of the soft function. We turn to this next. 

In the spirit of the form factor scheme that we employ here, we will associate a jet function in (|10[) for each of the 
hard partons, both massless and massive. The number of hard partons in the same underlying process determines 
also the number of jet functions in the decomposition of an observable ©. To derive the expressions for various jet 
functions in the following we use their process independence to either calculate them directly or to extract them from 
known results. 



D. Jet functions for incoming partons 



We start with the jet function Ji n for an initial state hard parton (quark or gluon) which is basic for all hadron 
collider processes. To this end, we can use the well known results from Drell-Yan vector boson or Higgs boson 
production. In these reactions, Eq. © takes the form 

<j p (N, Q) = [J?(N, Q)] 2 H(Q)S(N, Q) + 0(l/N) , (13) 

for P € (q «-> DY,^ Higgs). In the two processes, the hard scale Q is simply the virtuality of the outgoing color 
singlet vector boson in DY or the mass of the Higgs boson. Since in these two reactions exactly two hard colored 
partons are involved, the hard and soft functions are just lxl matrices, i.e. the color structure is trivial. Upon 
setting ln(/i 2 /s) = in the results of Ref. 113 1 , it follows that the soft-anomalous dimension matrix vanishes through 
two- loops (and possibly to all orders [l3L Il4. llH fl6l|). Therefore, in Drell-Yan and Higgs boson production we have 
simply S(N, Q) — 1. Thus, the jet function J; n is simply the square root of the corresponding Sudakov factors, see 
for example Refs. [22l. |23||: 

1 rl r N-l _ i f Al-x) 2 Q 2 , 2 \ 

\*J?{N,Q) = \ / dx X - / %2A P (a s [ q 2 ])+D P (a s [(l-xfQ 2 }) . (14) 

1 Jo l — x y j M 2 q j 

The functions Ap,P — (q,g) and Dp,P — ((?,<?) are currently known through three loops ([HI, Hf| and (22l. [26l. [27j 
respectively). The factorization scale fj,F appearing in Eq. (|14[) is related to the factorization of the non-perturbative 
parton distributions, assumed to be defined in the MS scheme. Utilizing a perturbative distribution function [28l. [29l. |30| 
one can also extend that result to processes initiated by massive partons [l8[. 

The derivation of the jet functions for final state hard partons is more involved since these depend on the definition 
of the observable. Similarly to Drcll-Yan/Higgs, one can use the vanishing of the soft anomalous dimension matrix 
(and thus the absence of non-trivial soft-gluon correlations) in any process involving two hard colored partons in order 
to extract various jet factors. For example, jet functions for "observed" outgoing hard partons (fragmentation) can 
be derived from semi-inclusive e + e~ annihilation to hadrons [3 it l32l . l33l . l34j |. Extension to the massive case can be 
done in a fashion similar to the case of Drell-Yan discussed above. 

Of particular interest to us in this work are observables with inclusive final states; a very well known example is 
inclusive DIS [U, H3| which can be treated similarly to Drell-Yan and e + e~, as discussed above. We are furthermore 
interested in processes with non-trivial color correlations, like the resummation of soft-gluons at NNLL in ti hadro- 
production. In order to calculate all jet factors that enter that observable we need to first specify the soft anomalous 
dimension matrix in this process which is done in section IIVI The calculation of the final state jet factors and the 
final result for the cross section are relegated to section [V] 

Finally we would like to comment on the process independence of the various jet factors discussed above. In 
principle, the presence of a process dependent hard scale Q indicates process dependence of the whole result. What is 
process independent is the functional form of the corresponding jet functions, while the dependence of the hard scale 
should be thought of as a sort of functional argument related to the phase-space for soft-gluon radiation available in 
the given process. Therefore in different processes the "argument" of the jet functions will in general be different but 
their functional form stays the same. 



III. FROM DIFFERENTIAL TO INCLUSIVE CROSS SECTIONS 



The resummed partonic hard scattering function lup(z) at fixed invariant mass is found from its moments with 
respect to z — M 2 / s, with M the pair invariant mass and s the partonic center of mass energy squared. The fully 
inclusive hard scattering cross section is then found by integrating over M, or equivalently, over z, and the result is 
a function of 



only. We must also integrate over the center-of-mass scattering angle (equivalently, fj above), but as we shall see, this 
does not affect our reasoning, and we suppress this integral for simplicity of notation. In expressing our results, we 
will find it useful to note that the ratio of pair and particle masses obeys the relation 
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Wo denote the inclusive cross section at parton level as <Jp(p, s). It is related to the differential cross section in z 
at fixed s by simply integrating over z, with lower limit p, 



a P (p,s) 



dz uop{z, M, to) . 



(17) 



We have observed that the singular dependence of ujp on z can be found in turn by an inverse transform, 

ujp{z, M, to) = / — Z - N lu]5 s (N, M, to) . (18) 
J 2m 

The z-resummed cross section is taken at fixed pair invariant mass M, and therefore fixed velocity in the hard- 
scattering cm., 



P. 



M 



1 - 



Am 1 
If 2 



= 1 - 



(19) 



Our goal is to relate the expression for tu res (N, M,m) at fixed M to the inclusive resummed cross section ap(p,s) 
with respect to p, as an inverse transform from moment space in terms of that variable. 

The resummed expression for u>(N) is given in Eq. For the following analysis, we make a slight change in 
notation, and recognize that the hard-scattering function H, which describes the short-distance part of the cross 
section, vanishes linearly in the center-of-mass velocity at absolute threshold, and lowest order in a s . Since (3m 
depends only on the ratio m/M, this quantity is fixed for Mellin moments with respect to z. It will, however, be 
important for ap. To make this trivial but important factor explicit, we change the notation slightly, and write 



H 



M 2 to 2 



p M 2 



so that Eq. ([6]) becomes 

^ cs (iV, M, m) = M 11 J, (AO Tr 



H 



p M 5 

Z /i z 



,&s(p 2 



~M 2 ~'m 2 ~' as[p 



(20) 



(21) 



where we have represented the jet functions schematically. We emphasize that moments in z at fixed M are equivalent 
to integrals over s. They thus leave (3m, but not p, fixed. Dependence on p enters only after the integral over M 2 , or 
equivalently z, as in Eq. (flT)) . 

We now relate the inclusive cross section, which depends on p, to the fixed-M 2 z moments of u> by 



a P (p, S ) = J dzJ^L 



z~ N oj^ s (N,M 2 ,m 2 ) 



H 



P M 2 2 , 
-, — ,a s (n , 

Z pi A 



(22) 



where in the second form we have observed that uip(z) vanishes for z > 1, and have exchanged the N and z integrals. 
We note that when we choose /i = cm, with c some constant of order unity, the full right-hand side of Eq. (|2"2"]) depends 
on M only through logarithms of the ratio 4m 2 /M 2 = p/z, see Eq. (fT6)) . Making such a choice of renormalization 
scale and changing variable from z to 



P 



we derive the desired form of an inverse transfrom, 



dN 



a P (p,s) = I ^-p- N+1 a P (N,m) 



with 



/>oo 

a P (N, to) = / r N Vl^F 1 II UN, e, a.) Tr [H p a s (p 2 )) S p (iV 2 £, a s )) 

1 



(23) 



(24) 



(25) 
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Here we have simplified the notation for the arguments of the jet, hard and soft functions somewhat to emphasize 
their ^-dependence. The scale of the running coupling is, as indicated above, of order of the quark mass, m. The 
relationship between the iV-dependence of the resummed cross section at fixed M in ujp s (N) and in <rp(N, m) can be 
found readily in the large- TV limit, by noting that the integral over £ in (f^5|) is dominated by the factor 

r N „ e NXn{l-f M ) ^ (26) 

which forces M 2 towards 4m 2 . Center-of-mass velocities (3m are thus exponentially suppressed. Corre- 

spondingly, the scale of energy evolution in the soft cross section, Eq. {7} is over an interval from m to m/N > mf3 2 M . 
For this range of energies, the evolution variable fi' in |[7J) is larger than the kinetic energy of the pair in their center 
of mass. We shall assume below that, as suggested in Ref. [35| . radiation in this energy range decouples from the pair, 
whose interactions give rise to Coulomb enhancements that appear as inverse powers of 0m- In the soft anomalous 
dimension matrix appropriate to this range of energies, the pair of heavy quark eikonals is effectively replaced by a 
singlet or octet eikonal line, with a separate term that describes the evolution of the pair. This approximation results 
in a smooth limit at absolute threshold (3m = 0. [6ll 

Corrections due to the logarithmic ^-dependence in the jet and soft factors are suppressed by inverse powers of 
N . Up to such corrections, the result is the Born cross section for heavy quark production in process P times the 
remaining jet, hard and soft functions, which we write as 

a P (N, m) = a£ oin (N) ]J J t (N, 1, a a ) Tr 

i 

where the hat on the hard matrix indicates that have factored out the A^-dependence of the Born cross section, which 
behaves at leading power in N as 

it^f^^vn^^ ( 1+o (£0 • (28) 

Eq. Ip7|) is the form that we will use below. As suggested above, we will evaluate the soft function S P (N, 1, a s ) using 
Eq. computed with a soft anomalous dimension appropriate to the range energy range m > [/ > p^rn. 



H p (l,a s ( M 2 ))S p (iV,l,a s ) 



1 + 



(27) 



IV. THE TWO-LOOP ANOMALOUS DIMENSION MATRIX AT ABSOLUTE THRESHOLD 



In this section we derive the relevant result for the two-loop anomalous dimension matrix T$- We also show that 
these results are enough to predict the full pole structure of the two-loop color averaged amplitudes for qq — ► ti. 

The one- loop massive anomalous dimension matrix for an am plit ude with n colored partons, J\f m of which are 
massive and with equal mass m has been known for some time (El. |38|: 



(i#)=l 



where = (j>i +Pj) 2 and cry = 2pi ■ pj 





( a ij 








K 


2pi ■ Pj 


— Sij - 


2 

- m' — 


2 




Xij 


Sij 


(1 


— Xij) 2 



Ti ■ Ti 



In ( 1 + xj 



2x 2 
1 - x 2 



ln(xi 



(29) 



mf — m 2 (with rriij — {0, m}). The space-like variables x^ read [39j |: 



1 - 



- 1 



1 - 



(30) 



~Vij 



where sa is 



when, in the unphysical space-like kinematics, all invariants Sy < 0. In specific applications some of them have to be 
continued to time-like kinematics. This can be done with the help of the replacement Xi 
now in the physical region Sy > 4m 2 and the 
such that they satisfy = an d can be 

for quarks or gluons. The index i labels the leg 
In parallel to the two-loop massless case [Ti 
from 2- and 3-eikonal (3E) contributions, i.e. 
three) external hard partons. 



•'time-like" variable is < < 1. The color generators Tj are defined 
either in the fundamental or adjoint representation of the color group 

where the generator is inserted; see also appendix [X] for more details. 
|, the two-loop massive anomalous dimension matrix is built up 

configurations where soft gluons are exchanged between two (resp. 



,(2) 



E 

(i#)=l 



Ti 



Ti-TjP i3 



(2) 



3E terms , 



(31) 



(i^j)&Mm 
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where, as for the massless case, K = (67/18 — it 2 /6)Ca — (5/9)Nl- Even before treating the 3E terms, we see that at 
two loops exchanges involving two eikonals take on the same color structure as in the one-loop anomalous dimension, 
Eq. ([291). 

Even without using explicit forms for the 3E contributions [36| , we have adequate information to study the behavior 
of the soft anomalous dimension in the range [i > f3 2 m, subject to our assumption of factorization, as discussed above. 
In [f| it was observed, for example, that 3E contributions to the reactions qq — > QQ and gg — > QQ vanish, either 
identically (when two eikonal lines are massless) or at u = t (for two massive eikonals). Given our assumption of the 
decoupling of soft radiation from the dynamics of the pair in the range /i' > /3 2 m, we may extend the anomalous 
dimension, appropriate for this range of energies to absolute threshold, (3 = y/l — Am 2 /s — > 0. To the order at 
which we work, power singularities in (3 associated with independent evolution of the pair of heavy eikonal lines 
cancel in the soft function S of Eq. ([7]). It is this simplification that enables us to present a full expression for the 
threshold-resummed inclusive cross section at NNLL. To derive this result, we need only the 2E diagrams. 

The 2E, dipole-type contributions, can be readily derived in complete generality. Here we note first that Eq. (j3"Tj) 
reproduces the known massless result for a massless dipole. Second, it reflects the fact that, similarly to the one loop 
case, the mixed corrections between massive and massless legs do not produce any power corrections in m 2 /sy. To 
verify that one-mass dipoles do not involve additional power corrections, we have repeated for this case the arguments 

(2) 

given below for the derivation of the function py ' . In that check we have used the recent two-loop calculation of the 
heavy-to-light form factor in QCD [10, [H|; see also [12, IH, Si|. For partial checks we have made use of the packages 
HPL 45J and FIESTA [46]. Note also that the absence of power corrections in the mass in the one-mass dipoles is 
related to the choice of the variable = 2pi ■ pj instead of = (pj + Pj) 2 in the first term of Eqs. (|29I31[) . 

(2) 

Finally, we explain how the functions py can be determined. The dependence on the indices of the function 

P^ in Eq. (|3ip is only through the corresponding kinematical invariant Sy, i.e. P^ = P^ 2 '{sij), and the dependence 
on Sij enters through the variable Xij defined in Eq. (|30p . That implies its functional form is universal and therefore can 
be extracted from the simplest two-loop amplitude with n = 2: the two-loop massive vector form factor F\ (7* - 
In this case Eq. @U) simplifies to 



r?W2) - -c F f 



(32) 



where a — s — 2m 2 and s = (pi + P2) 2 < 0. Of course, there are no 3E contributions for n = 2. 

As we remarked above, the soft anomalous dimension matrix is defined only up to a term proportional to the unit 
matrix. In this context, we can use that ambiguity in the definition of the soft function to define it through the 
condition H = 1 in the factorization of the form factor F\ = J ■ S ■ H following from Eq. (JTUJ) . From the known results 
for the form factor Ft and the jet function J [l|| and taking into account Eq. (fT2")) we derive: 

p(2) = K p {l) + p(2),n. ^ (33) 

where, similarly to the definition of P^ in Eq. |3T|) . the function equals the term in the square brackets in 
Eq. (f2"9"]) . The presence of the term p( 2 )' m in the above equation indicates that the property of the two loop massless 
amplitudes = K/2 is broken in the massive case by power corrections of the mass. The function p( 2 )> m 

is given by 

X2{1 + X2) ln 3 (x)+x 2 (l-x 2 )ln 2 (x) 



3 

+ (-(l-x 4 )\n(l-x 2 )+x 2 (l + x 2 )C2)Hx)+x 2 (l-x 2 )C2 + 2x 2 (Jj, (34) 

where is the Riemann zeta function: (2 = tt 2 /6, £3 = 1.202057 .... 

The function p( 2 )> m (and in particular its real part) does not vanish at threshold; for example, for a time- like 
argument x = — (1 — j3)/(l + 0) + is with [3 — y/T— AmP/s and Q < < 1, this limit is 

P^(x ~ -1 + is) = l ~^C A + g - £) l jc A + 0((3) . (35) 
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This result contains, as usual, a Coulomb enhancement in its imaginary part, which reflects the pair's internal 
evolution. 

Combining the results above, the two-loop soft anomalous dimension matrix for the two-to-two quark- and gluon- 
initiated reactions (see Eqs. (|A1[) and (|A2p ) takes the following form close to absolute threshold j3 — ► 0, 

4 2) = y r«+T3-T 4 P 3 ^ m , (for/?->0), (36) 

with pj$' m given by Eq. (|33|) . The explicit results for the matrices and T3 ■ T4 can be found in appendix lAl 

The mass-dependent soft anomalous dimension of Eq. (|32|) for processes with the color structure of the form factor 
was derived in Ref. 0] , using a slightly different scheme for the soft function. Eq. (f3"Tj) for soft matrices of arbitrary 
n-point amplitudes involving massive colored particles was presented in Ref . [8j . To determine the analogue of the 
function p( 2 )' m , the authors of that reference have utilized the results of Refs. (48| and We find agreement between 
our p( 2 )' m and (—1 times) the function appearing in Eq. (15) of version 2 of the arXiv preprint of Ref. [8]. The fact 
that the results of Ref. Q reproduce the IR poles of the massive form factor (which we have used to extract the 
function p( 2 )> m ) implies a non-trivial consistency between our results and the results of Refs. 

Moreover we have performed for the first time a truly non-trivial check of Eq. (|3 1 [) as a whole, by predicting the IR 
poles of the squared two- loop qq — > QQ amplitude and comparing them to the numerical calculation of Ref. [49f | . We 
have found a complete agreement between the predictions of our formalism and the color-averaged squared amplitude 
at two loops. 

In order to be able to make this prediction, we have noticed that in squaring the amplitude and summing/averaging 
over colors, any 3I?-type contributions in Eq. (f3Tj) with color structure f abc T a T b T c would vanish simply by color 
projection, and would not contribute to the squared amplitude at this level. Thus the calculation in question does 
not test for the presence of such terms at the amplitude level. 

The setup of our prediction is as follows: the amplitude M, multiplied by the Born amplitude and summed over 
spin/color, can be expanded in the coupling a s = a s {n 2 ) / {2tt) as 

M = AfW (e) + a s M {1) (e) + a 2 s M^ (e) + O {afj . (37) 

The factorization properties of on-shell amplitudes detailed in section Hi CI give the following prediction for the poles 
of the amplitude M through two loops 

M«(e) = + J«| M< >+O(e°) 

(38) 

The function J(e) represents the product of the four e-dependent jet functions corresponding to the two incoming 
massless and two outgoing massive fermions (see Refs. [50] for more details) and has an expansion in a s analogous to 
Eq. (|37[) . Similarly, Ti and T2 are the expansion in terms of the coupling a s of the anomalous dimension matrix T$ 
given through Eqs. (|29I31|) . 

Predicting all two-loop poles in the squared amplitude requires also the one-loop amplitude for the same process 
evaluated to sufficiently high order in e. We have calculated them separately. Our results will provide a non-trivial 
check on future extension of the results of Refs. |51| to the analytic calculation of the non-planar contributions in this 
process. 

V. NNLL RESUMMATION FOR TOTAL it HADRO-PRODUCTION CROSS SECTION 

We can summarize the results of sections [TT] and Mil for the resummed partonic total-inclusive cross section for ti 
pair production at hadron colliders in moment space by 

<j p (N,m 2 ,n 2 )=a p ma (N) [J P (N, m 2 , M 2 )] 2 [ J incl (JV, m 2 , M 2 )] 2 Tr [u p (m 2 , fi 2 )S p (N, m 2 , M 2 )] + 0(1/N) . (39) 

The index P = (q, g) labels the jet functions as well as the two reactions qq — > ti and gg — * tt. The factoriza- 
tion/renormalization scales are denoted by \i. The jet factors J p a are given in Eq. (fH|) with Q 2 = 4m 2 , where m is 




11 



the pole mass of the top quark. The only factors remaining to compute in (f39|) are functions J; nc i for the final-state 
jets and the one-loop correction to the soft function at TV = 1 in Eq. ([7]), which serves as a boundary condition for 
the evolved soft function. We turn first to the outgoing jets. 

The outgoing jet functions are specified by our choice of the form factor scheme, as described in Sec. Ill CI Their 
virtual contributions precisely cancel the terms subtracted from the soft anomalous dimension matrix T$ in Eq. 
(|29|) . and hence in the soft function S appearing in Eq. (|39|) . Specifically, we have subtracted those soft singularities 
corresponding to the low-mass limit of the outgoing legs. For a completely inclusive observable, like the total inclusive 
cross section, such factorization is not strictly necessary. The resulting expression, however, Eq. (|39|) . provides a 
unified threshold resummation for the total inclusive cross section and for the cross section at measured pair invariant 
mass s > 4m 2 , including the limit s 3> m 2 , where logarithms of the heavy quark mass can be important. 

As described above, the outgoing jet function can be constructed directly from the exponentiation of its infrared 
singularities in the low-mass limit, and therefore is of the form, 



Jind(N,m 2 ,fi 2 ) = expi ^ J dx— 



N-l 



- 1 



1 



r inc i (a s [V(l-i) 2 ]) 



(40) 



The jet anomalous dimension r; nc i is proportional to the unit matrix in color space. Specifically, it is given by the 
single poles of the logarithm of the massive quark form factor in the small mass limit, which defines the jet factor of a 
massive line in an amplitude [181 ]. and which we have adapted here for the form factor scheme. The explicit expression 
for ri nc i through two loops is 



a s {n 2 



C F 



-1 - In 



K 



-C F 



-1 - In 



Ca-l 



CfCa 



(41) 



The non-logarithmic part of r; nc i can be naturally expressed in terms of the anomalous dimensions G, K needed for 
the exponentiation of the massive form factor to NNLL [l8| , 



r inc i(non - log term) 



a s (/j, 2 ) 1 



Gi' 



(0) 



^ 2 ) 



1 

I 2 " 



G 



(o) 



K 2 



(42) 



The functions G and K are defined in [18[ as expansions of a s /(47r), hence the additional powers of 1/4 in the 
equation above, and with Nl active flavors. The function G is e-dependent: G n — J2i=o Gn'f 1 , and it equals (minus) 
the function G in the massless form factor [l^, [2(| [2l|. The origin of the term /3qG^ can be understood along the 
lines of Ref. [H. 

The last step remaining is the derivation of the boundary condition (1, Pi ■ (3j) for the soft function S, see section 
IIIBI We recall that the boundary condition is uniquely defined once the form factor scheme has been adopted. To 
extract it, we need to calculate the total inclusive cross section in the eikonal approximation. After the appropriate 
cikonal jet functions have been factored out (see Ref. [1]) we are left with the desired boundary condition. 

The required one loop calculation is in fact quite straightforward. To that end one can use the factorization in 
the soft limit of the squared one-gluon real emission amplitude into the square Sij of the eikonal current and the 
Born amplitude with appropriate insertions of the color operators Tj ■ Tj summed over all pairs of legs (i, j); see, for 
example, Ref. [53| for details. Combining matrix element factorization with the factorization of phase space in the soft 



limit, we arrive at a 



(1) ,real 



ikonai — J2i j=i Born y x Iij . We label the legs according to the momenta of the hard partons; 
see Eq. (jA2l) . The functions 1^ are simply the integrals of the eikonal current squared over the phase space of the 
emitted gluon. While the integrand is scaleless by construction, the integrals do not vanish because we integrate up 
to the maximal energy E max available to the emitted gluon in the partonic cm. system. Their expressions read: 



AE 2 

max 



J, 



17 > 



where : J,;. 



2 2 - 2e 7r 1 - 



E, 



e dE n / dn 



(Pi ■ g)(pj -9) 



(43) 



Working out the color algebra we get the following result for the one-loop real-emission contribution for the reaction 
99 ~ ¥ QQ (which covers the general case), 



(l),real 
eikonal 



0"Bo 



a a (n 2 ) 



AE 2 



-2 (C F (J34 — J 33 ) + Ca J 




Ca {J 12 + J34 — 2J13) 










(i 


1 












(44) 
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where cr Born above is a diagonal matrix; we work in the singlet-octet basis given in the appendix [5] 

Eq. (|44[) is derived in the back-to-back scattering configuration where u = t and holds for any (3 — yl — 4m 2 /s. 
Nicely, the result is diagonal and the two octets are degenerate. To complete this result one has to add the cor- 
responding virtual corrections. Since they are all scaleless and thus vanish in dimensional regularization, the only 
contributions comes from their counterterms. 

The result simplifies significantly if one takes the limit (3 — > which is relevant for the resummation of the total 
cross section. In that limit (i.e. working up to powers of (3) we can set (3 = everywhere in the integrals Jij. In this 
limit we have J13 = J12/2 and J33 = J34 as well as 4:E 2 lax = s/3 4 . The two independent integrals read: 



J12 
•hi 



e e7E 1 T(-e)r(l - e) 
r(l-e)2e r(l-2e) 

e^ E 1 2 2£ V^r(l - e] 
r(l-e)4e r(3/2-e) 



J_ 7T 2 

2^ + ~8 



O(e). 



(45) 
(46) 



After subtracting the eikonal jets in such a way that the singlet eigenvalue vanishes (i.e., to reproduce the well known 
result from Drell-Yan-type processes) we finally get, 



7 (1) 

eikonal 



C, 














f f y 


1 


1 


• 


v 4.E 2 / 

\ max / _ 




\ 








(47) 



The constant coefficient above follows from the constant term in Eq. ([46|) . We have verified that Eq. (j47|) correctly 
reproduces the ln/3 terms in the color-singlet color-octet difference of the total inclusive cross-section (see Eqs. (8-10) 
in Ref. H). 

Following the discussion in section [TTT1 we set the ratio fj, 2 /4E 2 lax to unity so that the logarithmic term in Eq. ([47]) 
vanishes. This way we get the following result for the boundary condition of the soft function in Eq. ([7| relevant for 
the resummation of the total heavy-pair cross-section, 



S(l,a s (Q 2 /7V 2 )) = S<°> 



= S<°> 



{Q 2 /N 2 



C A 










[I 


1 












C A 



7T 



(/i 2 ) /3 . (N 2 f 




(48) 



which as shown results in a term (3q (Ca/2) In TV in the NNLL result. The 0(a s ) correction appears only when the 
pair is produced in an octet configuration at short distances. The result for qq QQ follows from Eq. (|A16[) . 

We are now ready to combine our previous findings and present our result for the resummed heavy quark cross 
section in moment space up to NNLL. Working in the singlet /octet basis for the soft anomalous dimension T$, where 
it is diagonal [6j , the result for the resummed cross section for heavy pair hadroproduction reads 



a p (N, m 2 ,/i 2 ) 



Tr 



(Q 2 /N 2 



C A U 8 



exp < 


Jo 







,4m 2 (l-z) 2 ,2 




0{1/N,N 3 LL), (49) 



where Ilg projects onto the color octet states (see Eq. (|A16I) 'l . and where D p — contains single-logarithmic anomalous 

dimensions from both the jet and soft functions, in a color-diagonal form (see below). This expression is our central 
result. It may be cast in a more familiar form, by combining the constant piece of the soft function into the hard 
function, and generating the In N dependence in the soft function from a slightly modified version of the function 
D^—, which we denote as simply and which includes a new term proportional to a 2 /?o, 



QQ 



<7 P (N, TO 2 ,/i 2 ) 



= Tr 



H p (m 2 ,^i 2 ) exp 

im 2 (l-x) 2 ^ 2 




2A f 



C(l/iV,iV 3 LL), (50) 
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where 



Dj, = Dpl + 2Rer, 



1(^0^(11). (5!) 



The explicit form for D^— uses the results for Tp in section HVl and appendix [A"l In the reactions qq — > and 
.95 ~~ * QQ: which we label respectively by P £ (qq, 53), it reads through two loops: 

- ^(-^> ( s : ) + (^) 2 {*? ( ; : ) + (-< - ^ - -4) ( s : . ; 

' (52) 

where D q q = £>dy and D gg = -Dmggs [HI]. Corrections to Eq. ([50)1 begin, as indicated, at next-to-next-to-next-to 
leading logarithm, and are determined by the three loop soft anomalous dimension matrix and the inclusive soft 
function at two loops. 

The hard function H(m 2 , /j?) in Eq. (|50|) is known exactly through one loop [HH (see also Ref. [55[). Interestingly, 
the total cross section is the only tt observable for which at present the hard function is known beyond the leading 
order with full color dependence. We have used the degeneracy of the eigenvalues of the matrix Tg in the gluon 
fusion reaction (see also appendix [AJ to explicitly perform the trace over the degenerate octet ei gen value in Eq. (|50|) . 
Therefore the hard function H in this reaction is also a two-by-two matrix as computed in Ref. [54j | . 

The interplay of the jet anomalous dimension T- lnc \ with the soft matrix is quite interesting. As can be seen from 
the results in appendix \K\ through two loops r; nc i equals minus the singlet component of the anomalous dimension 
matrix Tg. Indeed, it is natural to expect that the anomalous dimension for producing a color singlet is the same 
as the one in Drell-Yan or Higgs production, a fact that has been anticipated already in Ref. Reproducing this 
property, without imposing it by hand, represents a very strong check on the consistency of our setup and results. 

The vanishing of the sum of the color singlet anomalous dimension and r; nc i is even more striking given the fact 
that they refer to very different kinematics: the former is related to the "very heavy mass" limit close to partonic 
threshold where the mass is comparable in value to the hard scale, while the former is extracted from the small mass 
limit where the mass is negligible compared to the hard scale. This result is quite intriguing, and it is clear that it 
is not accidental, as implied by the argument that the singlet anomalous dimension should not receive corrections 
beyond Drell-Yan/Higgs. Clearly, one can explore this property to relate the anomalous dimension in heavy flavor 
hadroproduction to the small-mass limit of the form factor beyond two loops assuming, of course, the findings of 
Ref. [y] extend to three loops. Combining the results in [l8[ and [47| one can cast the three-loop result in terms of 
only one unknown constant K3 (note that the function G3 is known from the massless form factor [2lL l56l|). 

The explicit form given in Eq. (|52|) is among the main results of this work. These anomalous dimensions provide 
explicit corrections to the form proposed in Ref. (57j . which was based on generalizing the proportionality between 
the one- and two- loop massive anomalous dimension matrices of the massless case [13| . The results of the present 
work as well as of Refs. [H, 0, HI provide the necessary corrections, which arise even for the special kinematics close to 
absolute threshold. 



VI. DISCUSSION AND OUTLOOK 



In the present paper we have extended the formalism for the resummation of soft gluon logarithms in cross sections 
with massive partons to the NNLL level. A central role in our construction is played by the massive two-loop soft 
anomalous dimension matrix for processes with n > 4 colored hard partons. In this paper we have presented the most 
general form of the so-called 2E (dipole) contributions. Combined with the results of Ref. @ this allows resummation 
in observables with special kinematics, like the total inclusive cross section for hadroproduction of a pair of heavy 
flavors. 

In our discussion we have detailed the relation between the soft function in an observable and the soft function in 
the corresponding factorized (virtual) amplitude. Following [f|, we have shown how the two are closely related for 
generic processes, and that the infrared poles of the amplitudes can be used to specify properties of the cross sections 
near partonic threshold, particularly in the form factor scheme defined and applied here. 

We have also explained how to construct the various jet factors needed for the completion of our threshold re- 
summation, and have used their process independence to derive initial-state partonic jets from the Drell-Yan vector 
production process. We have also used heuristic arguments to identify and construct natural jet factors needed for 
inclusive observables, like the total cross section for heavy pair production. Factorized jet functions, although not 
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strictly necessary, provide a form that can be extended to cross sections at measured pair mass, even far above absolute 
threshold. 

The most phenomenologically relevant application of our work is the total cross section for heavy pair production 
at hadron colliders. In this paper we have shown how to derive this quantity from resummed cross sections at fixed 
scattering kinematics. We have also given the exact result for the two-loop anomalous dimensions controlling the 
exponentiation of the NNLL terms in the cross section, and we have verified that even above threshold the part of 
the two-loop soft anomalous dimension constructed here is adequate to determine exactly the pole structure of the 
two-loop color-averaged amplitudes for top production through light quark annihilation. 

Our result provides not only the result for the resummed logs to NNLL but also a framework for studying the higher 
order effects in this observable and the associated theoretical uncertainties. We will provide a detailed numerical 
analysis in a dedicated publication. 

The formalism we have presented here can also be applied to more differential observables in heavy flavor production 
at hadron colliders, including cross sections at fixed rapidity and pair mass. To complete such studies one will, however, 
require explicit results, whether analytic or numerical, for the 3-E-type contributions to the anomalous dimension 
matrix involving two massive partons. 



APPENDIX A: COLOR BASES 



In this appendix we present calculations of the one- loop anomalous dimension in singlet-octet basis, and the eval- 
uation of the two-loop expression of Eq. (136)) . 



1. One-loop results in singlet-octet basis 

We work out the general result for the one-loop soft anomalous dimension matrices in the form factor scheme for 
the two reactions: 

«(Pi)?(Pa) -> Q{pz)Q_(Pi) (Al) 

g(pi)g(P2) Q(pMp±) (A2) 

where p\ = p\ = 0; p\= p\ = ra 2 . We define the invariants s — {p\ +P2) 2 , t\ = (pi — p^) 2 — m 2 = (j>2 — Pa) 2 — m 2 , u\ 
1 1 — p 4 ) 2 — m 2 = (j>2 — P3) 2 — m 2 . In the massive case 034 = 2p% ■ p^ — s — 2m 2 7^ s. 

We first consider the reaction (|Alj) . To evaluate the color matrices 1$ • Tj a color basis needs to be specified. The 
simplest one is: 

ci = <5i2<5 3 4 , c 2 = <5i 3 <5 24 , (A3) 

There are six combinations of the indices that need to be considered. Only three of them are different. Denoting 
these three color matrices by Ts, Tt and Tu defined through Ts = T\ ■ T2, etc., we get: 

m - - - /r (1) L (1) \ 

if 5 = 2T r T + 2U Tu + (So + S m + P) r s = 'ft J$ ' 

\ I 21 1 22 / 

where: 

rft = (2U-2T-S -S m -P)C F + (T-U)C A , 

rf 2 } = ±{2U-s -s m -p), 

= U-T, 

r$ = (2U -2T - So - S m - P)C F - (2U - So - S m - P)^ . (A4) 

The individual matrices Ts , Tt and Tu can be easily read off from the above equations. We have also included an 
overall minus sign in f s , fx as follows from the sign conventions for the color generators of Ref. [HI, |]58| : a generator 
Ti is multiplied by — 1 if it is inserted in a line that represents incoming quark or gluon or outgoing antiquark. 
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The expressions for the velocity factors U, T, So, S m and P follow directly from Eq. ([29)1 : 



U = \XL \ - — 

Ui 



T = ln , 5 = ln ^ +I7T, 



1S77 



In 



2m 2 



■in , P = P. 



(i) 



1-/9 
1 + 



(A5) 



We have labeled these functions according to their respective channel (s,t,u), and as to whether they refer to a 
massless or massive dipole. The function P collects all power corrections in the mass, i.e. in the massless case 
S m = Sq, P = 0. 

For physical applications to heavy flavor hadroproduction one chooses the s-channel singlet-octet color basis: 



1 



1 



vi 



Cl , v 8 



The change of basis for the anomalous dimension matrix in Eq. (|A4[) (or for any one of f s, f t, f u) is: 



-,(l),S.O. 



r/ ■ I t(1) t> 
xv • 1 § ■ xv. , 



where the transformation matrix R reads: 



R 



1 




1 

2 



R 



- 1 - N 

2 



The matrix elements of the matrix T 



r (i),S.O. _ , q 
1 11 — - {00 



(i),s.o. 



read: 



,(l),S.O. 



r 



^(l),S.O. 
- 21 

(l),S.O. 
22 



2(U-T), 

(4C7 - 4T - S - S m -P)C F - (4C7 - 2T - 5 - 5 ro - P) 



C A 



(A6) 



(A7) 



(A8) 



(A9) 



Rearranging the above expressions we arrive at the following expression for the anomalous dimension matrix rg' ,S '°' 
in the singlet-octet color basis: 

r (i),S.O 

r 
r 



11 

(l),S.O. 
12 

(l),S.O. 









c F 




--Or) 


— Lp — in 









21 



.(l).S-O. 
22 



hi 

2 In 
4 In 
2 In 



C F 



ui J Ca 

h 

Ul 

h 

Ul 

h 



2 In 
2 In 



_ui_ 

s 



In 

-In 



m ~ 

s 

m' 
s 



Lp ~ ITT 
' Lp — ITT 



C F 

Ca 
2 



We have kept the traditional notations and defined: 

l + P 2 



Lp = 



2(3 



In 



1 + /3 



(A10) 



(All) 



The result for the anomalous dimension matrix in rg' ,S '°' in Eq. (j A10|) agrees with the one derived first in Ref. Q 
provided we add to the above result (|A10j) the term (— ln(m 2 //Lt 2 ) — 1)C F 1 and set /i 2 = s. The addition of this (color 
diagonal) term corresponds to working in a scheme where the two jet factors for the heavy quark and antiquark are 
absorbed into the soft function (see the discussion in section [Vj . 



16 



Finally, we give the expression for the real part of the anomalous dimension matrix rW' S O - a t absolute threshold. 
To that end we set u% = tx = — s/2 as well as ln(^ 2 /s) = 0. Up to corrections ~ O{0) the result reads: 



ReT 



(l),S.O. 



= Ct, 



= -r 



(i) 



1 o 
o 1 



i 



Ca 
2 



1 
1 




1 



Ca 
2 




1 



(A12) 



2. Evaluation of Eq. flMJ 

The threshold limit of the two- loop matrix Tg in the basis (|A6|) can be easily derived from Eq. ([36"]) by using that 

Ts-n = tf° =(-^ ° V (A13) 



-C f , 2 



(2) 

Up to corrections ~ 0(/3), the result for reads: 

Reri?-° = -iffi ( J ° ) - [Sf- + 6^i0l) ( ° °) , (A14) 

where and lj^j are given in Eq. (|41) . 

The calculation of the gluon fusion reaction Eq. (|A2[) can be done in the same way as in the quark-antiquark 
annihilation one described above. The appropriate s-channel singlet-octet color basis is: 

vi - S 12 5 34 , vg = d 12c T 3 c 4 , v a s = if 12c TI A . (A15) 

fll S O (2) S O 

A direct calculation shows that in the limit (3 — > the matrices T3 • T4, Re ' and Re r^ ; ' ' ' for the gluon fusion 
reaction can be obtained from the corresponding matrices in the quark-antiquark initiated one Eqs. (|A13IA12IA14[) 
with the help of the simple replacements: 

I 2 x2^l3x3 and ^jj^j_>|010|. (A16) 

The results of Ref. [59] can be used to simplify the calculations. 
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